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Debuter sur cas où les deux corps sont en orbite circulaire, ces points représentent les endroits où un troisième corps de masse 
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EAAVIKà En mécanique céleste, il est un sujet qui a passionné de nombreux mathématiciens : c'est le problème dit des trois 
English corps. Newton, après avoir énoncé sa loi qui exprime que « les corps s'attirent avec une force proportionnelle au 
so produit de leur masse et inversement proportionnelle au carré de la distance de leurs centres », a cherché à décrire 
spaño | a | 
a le comportement de trois corps sans y parvenir. Il faut attendre le mathématicien Joseph-Louis Lagrange qui, en 
Euskara 1772, étudia le cas d'un petit corps, de masse négligeable (ce qu'on appelle aujourd'hui corps d'épreuve ou 
E particule-test), soumis à l'attraction de deux plus gros : le Soleil et, par exemple, une planète. Il découvrit qu'il existait 
Suomi des positions d'équilibre pour le petit corps, des endroits où toutes les forces se compensent. 
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aS] (échelle non respectée) : 
Lietuvių | | a | 
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Modifier les liens 


comparable à celle entre L4 et ce corps. 

e L, : sur la ligne définie par les deux masses, au-delà de la plus grande. Dans le cas où l'un des deux corps est 
notablement moins massif que l'autre, la distance entre L3 et le corps massif est comparable avec celle entre les 
deux corps. 


e L, et L; : sur les sommets des deux triangles équilatéraux 
dont la base est formée par les deux masses. Sans qu'il y ait 
de consensus précis, L4 est celui des deux points en avance 
sur l'orbite de la plus petite des masses, dans son orbite 
autour de la grande, et L5 est en retard. Ces points sont 
parfois appelés points de Lagrange triangulaires ou points 
Troyens, du fait que c'est le lieu où se trouvent les astéroïdes 
troyens du système Soleil-Jupiter. Contrairement aux trois 
premiers points, ces deux derniers ne dépendent pas des 
masses relatives des deux autres corps. 


Calcul de la position des points de 
Lag range [modifier le code] 


Les deuxderniers points de Lagrange forment =- 
avec les deux corps des triangles équilatéraux 


Le calcul de la position des points de Lagrange se fait en 
considérant l'équilibre d'un corps de masse négligeable entre le 
potentiel gravitationnel créé par deux corps en orbite et la force 
centrifuge. La position des points L4 et Ls peut être obtenue analytiquement. Celle des trois autres points L4 à L3 
s'obtient en résolvant numériquement ou éventuellement à l'aide d'un développement limité une équation algébrique. 
La position de ces trois points est donnée dans le tableau ci-dessous dans le cas où la masse d'un des deux corps 
(en l'occurrence le numéro 2) est négligeable devant l'autre, situé à une distance R du précédent. Les positions sont 
données le long de l'axe reliant les deux corps, dont l'origine est identifiée au centre de gravité du système, et dont 


l'orientation va du corps 1 au corps 2. Les quantités r et q dénotent respectivement la position du corps 2 sur l'axe et 


le rapport de la masse du corps le plus léger à la masse totale des deux corps. Enfin, on utilise la quantité £ définie 
par €= (q/ 3)18 


Point Position par rapport au centre de gravité du système 


L a+ R (-. a z de Le _ ot) 


9 
L ra+ R (e — = — Pf — oc) 
5 1127 
aaa S 3 Olh) 
Ls à 129 t 207367 + (a ) 


Dans la littérature, on trouve parfois des expressions quelque peu différentes, du fait que l'origine de l'axe est prise 
ailleurs que sur le centre de gravité, et que l'on utilise comme terme à la base du développement limité le rapport 


entre les deux masses plutôt que le rapport de la plus petite à la masse totale, c'est-à-dire que l'on utilise parfois la 
CEPE rs M q 
quantité q' définie par g = — = 7 ~~. 
M l-g 
Détail du calcul — Introduction 
Préliminaires 


On note M4 et m la masse des deux corps, la masse du premier étant supposée supérieure ou égale à celle du second. Les 
deux corps sont supposés être en orbite circulaire, leur séparation étant R Les deux corps orbitent autour de leur centre de 
gravité commun. On note 4 et r les distances algébriques des deux corps par rapport à un axe orienté du corps 1 au corps 2 
(c'est-à-dire que r va être négatif et > positif). Le centre de gravité est défini par l'équation 

M: Ti + Mao = D; 
avec par définition, de la distance R, 

Ta — Y] = R- 


Ces deux équations ont pour solution 


za ip 
LE M 
Wii 
EE R 
° M 


où on a noté M= M, + m la masse totale du système. 
Les deux corps orbitent l'un autour de l'autre à une vitesse angulaire w, dont la valeur est donnée par la troisième loi de Kepler : 
> GM 

o E 


G étant la constante de gravitation. 


LL 


Si l'on se place dans le référentiel tournant avec les deux corps, c'est-à-dire à la vitesse angulaire w, un corps immobile sera 
soumis, outre aux forces gravitationnelles des deux corps, à la force centrifuge. Si on note r le rayon vecteur de ce corps, la 
force centrifuge par unité de masse f à laquelle il sera soumis s'écrit 


http://frwikipedia.org/wiki/Point_de_Lagrange 


2/14 


Point de Lagrange — Wikipédia 


09/08/2014 


fe=ru 


Équation fondamentale 


La définition d'un point de Lagrange est que la somme des forces gravitationnelles et inertielles s'annule en ces points. En 
notant r le rayon vecteur du ou des points en question, on a ainsi 
GMi(r—-r1) Gmalr — rə) 2 
a le e D: 
|r — rall |r — rall 
les doubles barres indiquant que l'on prend la norme des vecteurs considérés. On remplace ensuite la vitesse angulaire w par 
sa valeur issue de la troisième loi de Kepler, ce qui donne 


GMi(r-r1) Gmar-r:) GMr 


Z a a — (} 
[r — ril? |r — ral? R? 
que l'on simplifie immédiatement par la constante de gravitation 
r— ri r — ro r 
M à NE D 
Trent Tr 


C'est la résolution de cette équation qui donne les différents points de Lagrange. 
Les deux cas à considérer 


La projection de cette équation perpendiculairement au plan de l'orbite, dont la normale est donnée par un vecteur noté p 
donne immédiatement 


r-2z = (0 
ce qui implique que l'ensemble des points de Lagrange est situé dans le plan de l'orbite. La résolution de l'équation se fait donc 
dans le plan orbital. Deux cas sont à considérer : 


e celui où l'on cherche un point le long de l'axe formé par les deux corps, 
e celui où l'on cherche un point en dehors de cet axe. 
Le second cas s'avère être le plus simple à étudier. 


Détail du calcul — Les points L; et Ls 
Cas des points L, et L; 


On suppose que le rayon vecteur r n'est pas parallèle à l'axe passant par les deux corps. On projette donc l'équation 
fondamentale perpendiculairement à cet axe, direction que l'on suppose définie par un vecteur noté Y Par définition, cette 


direction étant perpendiculaire à l'axe reliant les deux corps, on a 
Y- -rı= y- r50. 


L'équation fondamentale se réécrit donc 


ÿ y y y 1 T 

|r — r| [r — ral R 

Les termes en y . q Se simplifient, ce qui donne 
E 

On définit maintenant la direction g comme la perpendiculaire à r. Comme r n'est pas colinéaire à r4 et r,, les quantités 
FE T4 gne sont pas nulles. En projetant l'équation fondamentale le long de s, on obtient 

8 - ri CE Ta 

Mi = + Mo—— = 0 
lr — rall |r — rall 


Or, d'après le théorème de Thalès, les projections de r4 et r le long de 8 sont dans le même rapport que les projections de 
ces vecteurs le long de l'axe reliant les deux corps. Il s'ensuit que l'équation précédente peut se réécrire 


M È 0. 


“Ilr = ral 


Le barycentre des deux corps implique, comme wu précédemment, que 


[r-re 7 


Mi Ti + Mefa = D: 
La combinaison de cette équation et celle qui précède implique donc que les deux distances | | r— fr] | | et | | T— To | | sont 
identiques, leur valeur étant notée R : 

FAI elle 
En injectant ce résultat sur la projection le long de r, il vient alors 


M: To M 
"E PTE) 


En multipliant le tout par RÈ et en se souvenant que M est la somme des deux masses, on obtient finalement 


ce qui donne au final 
[r = rll = [lr — ra|| = E 


c'est-à-dire que les points cherchés forment un triangle équilatéral avec les deux corps du système. Ces triangles sont de plus 
inclus dans le plan orbital, ce qui donne deux points possibles, notés comme annoncé L, et Ls, étant situé de part et d'autre de 
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l'axe reliant les deux corps. 
En utilisant le théorème de Pythagore, la distance D de ces deux points de Lagrange du centre de gravité du système s'écrit 
2 2 2 
D? = (ri + R/2) + (V3R/2)> 
ce qui donne 


l 3 
Drd ru — Riral + a 


ce qui donne 

D? = Ë +ri — Rijril 
En utilisant le fait que fẹ — [ral + ro il vent 

2 2 2 

D'=r + |r| R = ri + raf 
La distance est donc supérieure aux distances de chacun des deux corps au centre de gravité du système. Ces points de 
Lagrange sont donc au-delà de l'orbite du corps le moins massif et ne sont pas strictement situés sur celle-ci, quoique ce soit 
quasiment le cas dans la limite où la masse du corps le plus léger devient négligeable par rapport à celle de son compagnon. 

Détail du calcul — Les points L4 à L, 

Cas des points L; à L, 
Dans le cas où l'on considère des points de Lagrange situés sur l'axe reliant les deux corps, trois sous-cas sont à considérer : 


1. Le cas où le ou les points sont entre les corps 1 et 2 ; 
2. Le cas où le ou les points sont à l'opposé du corps 2 par rapport au corps 1 ; 
3. Le cas où le ou les points sont à l'opposé du corps 1 par rapport au corps 2. 


Dans ces trois cas, l'équation fondamentale se réécrit de la façon suivante : 


e Cas 1 : la projection des deux forces sur l'axe a des signes opposés (la projection de la force exercée par le corps 1 est 
négative, celle de la force exercée par le corps 2 est positive), ce qui donne 


se = 
(r _ m (r = ra)“ Rè? 
avec 
FL TT 
e Cas 2: la projection des deux forces sur l'axe est négative, ce qui donne 
Mi 0 
(r = r1)° (r n ra)? R? 
avec 
TiS TY 
e Cas 3: la projection des deux forces sur l'axe est positive, ce qui donne 
E E E E) 
(r-r)? (r-r)? R? 
avec 
r < ry To 


Chacune de ces trois équations peut se ramener à une équation polynomiale du cinquième degré, pour laquelle il n'existe pas 
de solution analytique exacte, sauf cas particulier (comme celui des deux masses identiques par exemple). 


L'unicité des solutions dans chacun des trois cas se déduit du fait que l'équation à résoudre sur l'équilibre des forces dérive 
d'un potentiel U, donné par 


M; Tla 1 Mr? 
r= rill [rl 2 À 


Ce potentiel représente des pôles en 4 et », et correspond en dehors de ces valeurs à la somme de trois termes concawes et 
est donc localement concawe. Il ne possède donc qu'un seul extrême local dans chacun des domaines où il est défini, c'est-à- 
dire dans chacun des trois cas cités plus haut. 


U = 


Solutions pour L; à L} dans le cas où le rapport entre les masses est faible 
Forme réduite et solution dans le cas où le rapport entre les masses est faible 


Quand le rapport entre m et M} (ou entre m et M) est faible, on peut trouver une solution approchée pour la position de chacun 
des points en effectuant un développement limité à partir d'une solution approchée facile à trouver. Pour simplifier les notations, 
on effectue un changement d'échelle afin d'exprimer toutes les longueurs en unité de la séparation R et les masse en unité de 
la masse totale M On pose ainsi 


u=r/h 

et 
1,2 = r1,2/R 

et on définit le petit paramètre q par 
q = ma M 

à partir de quoi on peut exprimer 
M /M=1-q 


t = : 
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u =l- g 
Dans ce cas là, les trois équations écrites ci-dessus prennent la forme plus simple 
e Cas 1: 
Lei q 


D + a j) 
u+ uit 


avec 
-q< u<l-gq 
e Cas 2: la projection des deux forces sur l'axe est négative, ce qui donne 


l—g q 


(u+? (u-1+q)? ns 
avec 
1-g<u 
e Cas 3 i la projection des deux forces sur l'axe est positive, ce qui donne 
— 4 q 
CEE TES EN RSR 
avec 
u < —Ù 
Le point L, 


Quand la masse du corps 2 est négligeable, l'attraction de celui-ci est négligeable sauf si la particule d'épreuwe est très proche. 
Or, quand l'attraction du corps 2 est négligeable, l'équilibre entre l'attraction du corps 1 et la force centrifuge est tel que la 
distance du point d'équilibre est de l'ordre de R. Quand le point d'équilibre est situé à l'opposé du corps 2, on est dans le cas 
du point de Lagrange L3, qui est donc, en gros, situé à l'opposé du corps 2 par rapport au corps 1. Dans le cas contraire, on va 
donc supposer que le point d'équilibre est plutôt proche du corps 2 (et donc à nouveau situé à la distance R du corps 1), mais 
néanmoins suffisamment éloigné pour que l'attraction du corps 2 exercée sur la particule d'épreuve reste petite par rapport à 
celle du corps 1. On pose donc à partir de la forme réduite 


u=mm+te=l-g+e 


où ici £' est une quantité petite et négative (on suppose ici que le point est entre les deux corps). L'équation réduite se 
transforme alors en 


l-q 4 i 
e a LE = 0. 


On effectue un développement limité au premier ordre de l'attraction produite par le corps 1 : 
F q F 
-(-g)(-2e)+-;+1-g+e = ( 
Les termes en 1 - q se simplifient, et il reste 
1. 4 


Toujours en ne gardant que les termes d'ordre le plus bas en q, il vent 


On peut par la suite continuer le calcul, en développement l'écart du point au corps 2 en puissances de £". On pose ainsi 


u=1-g+e +r” + Oly 

L'équation fondamentale réduite donne alors 
l -g q Í F 
2 

—> ——— © + —— © Tl—-gtre +ze —=Ù 
On peut factoriser le second terme avec q / E2, que l'on peut remplacer par sa valeur, soit -3 &”. On obtient alors 

=(1-g){1+é+z) *-3e(1+xe) *+1-q+eé Ha =0 
On effectue ensuite un développement limité des deux premiers termes, au second ordre pour le premier et au premier ordre 
pour le suivant, ce qui donne 

(1-9) 26 — 25e + 37) — 3e (1- 9x) +1-q+e tre” = 0 
d'où on déduit que x vaut un tiers, ce qui donne 


1 m 
u=l-g+e + 3e +O(e) 


Le développement peut ensuite être continué suivant la même procédure. À l'ordre suivant, on a ainsi 


1 l3 
u = | — q + € + EVE = ~ + Oleh) 


Le point L, 
Le cas du point L, se résout exactement comme dans la section précédente, si ce n'est que le signe du second terme de 
l'équation fondamentale est négatif. On pose donc 

u=l-q+e 


€ étant cette fois-ci supposé petit et positif, et on a ainsi 
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l-q q4 
Ecoe a Li ni: 


La résolution à l'ordre le plus bas donne 
-(1 -g)(1- 2)- Z +1-qg+e=0 


qui après annulation des termes donne 
ka 
g2 


soit 


i i 
TXT à 
e= (:) =€ 


Cela correspond au signe près au même résultat que précédemment. La suite du développement de la solution se fait de 
même que précédemment. On part de 


u= 1— q4 e+ re 
et on injecte ce résultat dans l'équation fondamentale 
l—g q 
2 
— ————————— Hl Heare = 0Ù 
(l+e+re)? (e+re) j 


Comme précédemment, on transforme cette expression selon 


(1 = g)(1 — 2e — 20e + 3) — 3e(1 — 2re) +1- q+ e+ r =D 
ce que l'on résout en 


1 
T= > 


3 


soit 


1 | 
u=l-g+et+sé +O(e) 


Cette expression est identique à celle du premier point de Lagrange en remplaçant £' par €, mais ces deux points sont 
dissymétriques : comme le signe de £, £” change entre le point L4 et le point L>, la correction du second ordre, toujours positive, 
rapproche le point L du corps 2 alors qu'elle éloigne le point L> : les deux points ne sont plus à égale distance du corps 2. 
Pour la Terre, le rapport de masse est de 1 / 300 000, et £ est de l'ordre 0,01, ce qui place les deux points par rapport à la 
Terre à une distance d'environ un centième de la distance Terre-Soleil, soit dans les 1 500 000 kilomètres. Le terme de second 
ordre est de l'ordre d'un trente-millième de la distance Terre-Soleil, soit dans les 5 000 km. Le point L4 est donc environ 

10 000 km plus près de la Terre que ne l'est L. 


Enfin, on peut poursuivre le développement à l'ordre supérieur, ce qui donne, tous calculs faits 
1 1 
u=l-gte+ ze -2é + Of) 
3 g 
Le point L, 
Dans le cas 3, qui va correspondre au point L3, l'équation fondamentale s'écrit 
l-g q 
utg (u= lpg 
Comme le point est supposé au-delà du corps 1 par rapport au corps 2, il est plus proche du corps le plus massif, dont 


l'attraction va être prépondérante par rapport à l'autre corps. Dans la situation où l'on se place, le point recherché a donc sa 
position approximée par 


La solution approchée de cette équation est bien sûr 


u ~ —1 
Pour trouver les écarts à cette valeur, on écrit dans l'équation fondamentale 


u = —1 +2 
et on résout l'équation en prenant en compte les premiers termes en g. On obtient ainsi 
L= q 


M + — lHa 
(—1 +v +q) ü (1—1 +v+q)? Pi 


Les quantités v et q étant petites devant R, le premier terme s'écrit 


l—g 2 
"© = v—g)? (1 g)+ Av + g) 
TE +4) (1 — a) q) (1— g) + 2(v + 9) 
Le second teme étant négligeable par rapport au précédent (il est proportionnel à q), il peut s'approximer en 
q q 


(—-1—-1+v+q} TA 


En combinant l'ensemble de ces termes, on obtient 


1-q+2(v+0)+F-1+v=0 
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ce qui donne 

9 

-q + Su =Ü 
c'est-à-dire 


Oo öm ms 
name (3) 


On peut sans difficulté continuer ce calcul en posant désormais 


u = —=1 9 
u = —1 — l 21 + 
w étant cette fois proportionnel à g. L'équation fondamentale devient alors 
l-q q D 
À 
(—-1-Sg+w+g) (-1-1-£q+w+a) 
c'est-à-dire 
l-q q 5 
E R a D 
(1-4-0) (2-5-0) 
En développant cette expression au second ordre en g, on trouve 
w = 0+0O(g°) 
c'est-à-dire que west au plus en P. En refaisant le calcul dans ce cadre là, on trouwe finalement 
1127 , i 
H= 0 : 
w = ag736% TO) 


Il est rarement utile de pousser le calcul jusque là : dans une configuration Soleil-Planète, le demier terme correctif est au 
mieux de l'ordre 10°, puisque le plus gros rapport de masse Planète-Soleil, dans le cas de Jupiter est de l'ordre d'un millième. 
Le terme P est donc, pour Jupiter, de l'ordre d'un milliardième, ce qui, au w de la taille de son orbite, correspond à une 
correction d'une cinquantaine de mètres étant donné que la fraction en facteur de q° est de l'ordre d'un vingtième. Pour le 
système Terre-Soleil (distance de l'ordre de 150 millions de kilomètres, rapport de masse de l'ordre de 1 / 300 000), la dernière 
correction est une fraction de micron ! 


Sta b i | ité [modifier le code] 


Le calcul ci-dessus n'indique en rien si les points de Lagrange 
sont stables. La stabilité ou non de ces points est du reste peu 
intuitive. Dans le référentiel tournant avec les deux corps, une 
particule d'épreuve peut être vue comme soumise à un potentiel 
incluant la contribution gravitationnelle et celle de la force 
centrifuge. Ce potentiel, noté Q, s'écrit ainsi 


GM Gmo 1 


Tous les termes de ce potentiel sont négatifs et décroissent à 
mesure que l'on s'éloigne des masses (pour les deux premiers 
termes) ou du centre de gravité du système (pour le troisième). 
On peut ainsi montrer que les points de Lagrange L; et L5 sont 
des maxima locaux du potentiel Q (voir ci-dessous) et que les 
trois autres points sont des points selles. D'ordinaire, une 
position d'équilibre (déterminée par l'annulation des dérivées du 
potentiel) est stable uniquement si on se situe dans des minima 
locaux du potentiel. Cependant, étant donné que l'on est dans un 
référentiel tournant, le référentiel est non inertiel. Un objet se 
déplaçant dans ce référentiel, par exemple au voisinage d'une 
position d'équilibre, va être soumis à la force de Coriolis, et son 
mouvement ne dépend pas uniquement de la forme du potentiel. 
Pour étudier la stabilité des points de Lagrange, il faut donc tenir 


compte de la force de Coriolis. 


Pour calculer la stabilité des points de Lagrange, il faut ainsi Les points L4 et L5 bien que situés à des maxima 4 
étudier l'équation du mouvement d'un objet situé au voisinage du potentiel sont paradoxalement stables. L1, L2 et 
d'un de ces points. En notant R le vecteur de coordonnées ðX L3 qui sont des points-col sont instables. 
et òY donnant l'écart d'un tel objet à un des points de Lagrange 


(que l'on suppose confiné au plan orbital, l'équation du mouvement s'écrit 
åR = —2w AR +åf 
où ôf représente la force par unité de masse exercée sur l'objet. Cette force est petite du fait qu'au point de 


Lagrange la force (constituée d'une composante gravitationnelle et de la force centrifuge) est nulle et que l'on se 
place à proximité d'un tel point. Cette force peut se calculer en termes d'un développement limité. Par exemple, pour 
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la composante X, on a 
dfx dfx 


fx = fx(0) + OX 3Y 


Le premier terme correspond à la force s'exerçant au point de Lagrange, force qui est nulle par construction. Par 


ÒX + DY 


ailleurs, la force dérivant d'un potentiel, on peut exprimer les dérivées de la force en termes de dérivées secondes du 


potentiel : 
Po Pn 
IX = DTA — xoyr 
On peut ainsi exprimer l'équation du mouvement en termes des composantes selon 
: . 0 p’ 
OX = 2w0Y — IX — DC 
SY = —2wôX — 5X — ——5Y. 


OXOY OY? 


Ce groupe d'équation peut être mis sous la forme d'un système de quatre équations différentielles du premier ordre : 


OX 0 0 1 0 ÒX 
ð òY 0 0 0 1 Mid 
t| X | | -Ra Qay 0 -2 5X | 
5Y -Ray yy 2w 0 5Y 
où les dérivées partielles du potentiel Q ont été notées en indice précédé d'une virgule (par exemple, Q yx correspond 
à Q/0x) 


La stabilité du point de Lagrange considéré est obtenue en recherchant les solutions de cette équation. Pour cela, il 


plt. On va ainsi procéder à la diagonalisation de la matrice 


ci-dessus, que l'on notera A. Les valeurs propres trouvées vont correspondre aux quantités l ci-dessus, les écarts à 


suffit de trouver des solutions de type exponentielle, en 


la position d'équilibre étant alors une certaine combinaison d'au plus quatre exponentielles. La stabilité du système 
est assurée par le fait que les exponentielles ne croissent pas au cours du temps, c'est-à-dire que les quantités I 
sont soit négatives, soit complexes à parties réelles négatives. En fait, il n'est pas nécessaire de diagonaliser 
complètement la matrice, il suffit d'en trouver les valeurs propres, c'est-à-dire les solutions de l'équation 


det(A — AI) = 0 
Ce déterminant s'écrit 


— À 0 1 0 


0 — À 0 1 
=p =i À —2w | 
-Qey yy 2w —À 
et il vaut 
Y -+ um ZE vu FT dw A + (À re 2 yy = D 


Cette équation peut se ramener à une équation polynomiale du second ordre en A2. Les solutions de l'équation de 
départ sont donc deux couples de nombres opposés deux à deux. Pour que deux nombres opposés soient négatifs 
ou nuls ou alors de partie réelle négative ou nulle, il faut obligatoirement qu'ils soient des nombres imaginaires purs, 
donc que les solutions de l'équation en A2 soient réelles et négatives. Pour que ces solutions soient réelles, il faut 
donc que le discriminant soit positif, soit ici 


(Qas + on + 40°) — AN un — n°.) > 0 


Une fois ceci obtenu, il faut que les deux solutions réelles soient négatives, ce qui implique que simultanément leur 
somme soit négative et leur produit positif, ce qui implique 


Jm ne 2 + Au? > D 
one Ts > 


La stabilité d'un point de Lagrange est soumise à la réalisation de ces trois contraintes. Parmi ces contraintes, la 
dernière a une interprétation simple : le signe de la quantité {} (2 w y détermine si la position considérée 


est un extrémum local ou un point selle. En l'occurrence, la positivité de cette quantité implique qu'elle doive être un 
extrémum local, condition nécessaire mais non suffisante à la stabilité du point de Lagrange. Quand cette quantité est 
négative, on a un point selle et le point de Lagrange est instable. Par contre, de façon plus surprenante, un point de 
Lagrange peut être stable s'il correspond à un maximum local du potentiel, c'est-à-dire que Q xx + Q „y peut être 
négatif, pourvu que cette quantité ne dépasse pas la valeur critique de -4 w?. En pratique, c'est ce qui se produit 
dans certains cas pour les points de Lagrange L; et Ls. L'interprétation physique de cette situation est que la stabilité 
est alors assurée par la force de Coriolis. Un objet légèrement décalé d'un tel point va s'en éloigner dans un premier 
temps en façon radiale, avant de voir sa trajectoire incurvée par la force de Coriolis. Si le potentiel est partout 
décroissant autour du point, alors il est possible que la force de Coriolis force l'objet à tourner autour du point de 
Lagrange, à l'instar des nuages dans une dépression qui ne se dirigent pas vers le cœur de la dépression, mais sont 
contraints à une trajectoire circulaire autour de celui-ci. 
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Suite du calcul 

Préliminaire 
Pour étudier la stabilité des points de Lagrange, on va être amené à calculer les dérivées successives du potentiel. Ce potentiel 
fait intervenir la distance |r - r|. Il faut donc connaître les dérivées des différentes puissances d'une telle quantité. En 
coordonnées cartésiennes, cette quantité s'écrit 

rm Ve mr mea 
Sa dérivée par rapport à l'une des coordonnées x, y, z, collectivement notées x; s'écrit donc 
i 2( T; — Tii) Ti — Li; 
ð; | r— ri | = 9 z 5 5 = : 

(z= t1} + (y-n)? + (z-z)? |r- rl 

La dérivée d'une puissance p quelconque de cette grandeur est donc 

: —? 

a(r — ri") = piz: — mir — r?” 
En adaptant ce résultat aux dérivées secondes des quantités intervenant dans le potentiel, on a 

à. 1 -0y ain) mi) 

a Pen Lee 
rT ri rT Fi rT Fi 


ce qui, pour le potentiel complet, donne 


r GAA (z; — T1 i) (£; — T] i)GMi Gma (T; — La i) (£; — To ;)G M 3 
dj = TT RE o a a EEE —w di 
|r — ril |r — rl |r — ral |r — ral 
où ò; représente le symbole de Kronecker. C'est la valeur de ces dérivées partielles qu'il faut calculer pour déterminer la 
stabilité des différents points de Lagrange. C'est pour les points de Lagrange L; et L5 que ce calcul est le plus simple. 
Cas des points de Lagrange L; et L; 
Ces points se caractérisent par le fait que leur distance aux deux corps est identique et égale à R : 
|e — r| = |r — ra| = R 
Par ailleurs, on peut utiliser la troisième loi de Kepler pour passer des quantités du type G M/ R à w, et on connaît les 
coordonnées exactes des points de Lagrange. En évaluant les dérivées du potentiel aux points de Lagrange L4 ou Ls, on a 
1 V3 
Ti — Eli = hf, + — R 
2 2 
et 
1 V3 
Ti — Lo; = —- t, Ł+— R 
2 2 
le signe + s'appliquant pour L; et le signe - pour L4. Au final, la matrice recherchée a pour composantes 
3 3v3 
dN — y? TP. a5 4 S 2q) 
Fa (1 — 2g) T4 
Le déterminant de cette matrice vaut 
| 2T 2T 
; 2 — 23,2 B: 
det(3N) = DerDey — Voy = Te (1- (1 = 24)°) 2? = ali — 9) 
qui est toujours positif puisque q est confiné entre 0 et 1. Cette première condition de stabilité est établie. La seconde condition 
de stabilité s'écrit 
3 9 
FAIR (— = 7+4) hu 
| | 4 4 
quantité là encore positive. Enfin, le discriminant donne 
2 
2 2 2 2 
(Qas + y +407) A4 (Raray — a) = 2 (1 — 27q(1 — q)) = w? (27—27) 
La stabilité des deux points est, au final, déterminée par la positivité de la quantité 27q — 27q +1 Les zéros qa qp de 
ce polynôme sont donnés par la formule usuelle, qui ici indique 
| n 1 23 
fat — a Lilo 
z 2 2y 
Ce polynôme a ainsi des valeurs négatives sur la plage 
qE 1 Z 1 a 1 J 1 23 ~]0,03852, 0,96148]: Ainsi, la stabilité de ces deux points de Lagrange 
2 227 2 2727 
n'est assurée que si la plus petite masse n'excède pas 3,852 % de la masse totale, ou, de façon équivalente, que le rapport 
des deux masses n'excède pas 4,006 %*. 
Cette condition est vérifiée pour toutes les configurations de type Soleil-Planète (où q n'excède pas environ un millième pour 
Jupiter), ou pour le système Terre-Lune (où q est de l'ordre de 1/80, soit 1,25 %). 
Cas des points de Lagrange L; à L, 
Les trois points de Lagrange L4 à L; sont situés sur l'axe reliant les deux corps. Dans la formule qui donne les dérivées 
secondes, les quantités y; - y4 ; sont nulles, alors que leurs analogues en x s'identifient aux distances entre un des corps et le 
point de Lagrange considéré. En conséquence, la matrice des dérivées secondes s'écrit 
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2e — 2 Er — 0 
0 


à. () — TT] | r—ra| 
1] CM Gma PE 
Ir-r?  fr-r2ff 


Le terme Q xx est manifestement négatif. Le signe du déterminant de la matrice est déterminé par celui de Q „, : si ce demier 
est positif, alors le point de Lagrange est un point selle et il est instable. On peut réécrire ce terme en utilisant la troisième loi 


de Kepler : 
R? R? 
Ga =w = =t m 
A nE) 
Le cas de L, 


Le point de Lagrange L est situé entre les deux corps. Sa distance à ceux-ci, |r - r| et |r - r| est donc, à chaque fois, 
strictement inférieure à R. On a ainsi 
R? R? 


Q aiai (1 g)— © + g— © a egee 
ay lEs ( PTE |r — ra|? (( ) ) 
Cette quantité est donc strictement positive, ce qui assure que le déterminant est négatif, c'est-à-dire que L4 est un point selle, 
ce qui en fait un point instable. 

Le cas de L, et L; 


On pose, pour simplifier les notations, 


|r — r:l 
y = —— 
[rl 
SE 
On s'intéresse donc au signe de la quantité 
(2 1 1 
UU 1 
wa a ne 
W? ( be už 
soit 


Q 1 1 
“w (1-ġ|=-1 Zal 
se = (1-a) (+ }+4(& | 


sachant que u4 et u> sont reliés entre eux par le fait que leur différence est égale à 1 et qu'ils définissent un point de Lagrange, 
soit la relation 


1 1 |r 
—[ 1 — qì — — q7— — =. 
La distance du point de Lagrange au centre de gravité du système peut s'écrire, pour le point L», 
[r 
— =i +l- = W = g 
R (1 — q) q 
relations que l'on peut combiner en 
[r 
H Ta + (1 — gjur 


La position du point L> est donc donnée par 


On pose alors 
il 
A — 1 T To h e 7 
( q) F u) 


1 
tali 


On a donc, d'une part 


A+B=Ù0 

et, d'autre part 
(2 A B 
EE LESI Uo 


Autrement dit, 


AL ut 4 


Or, pour le point L>, on est situé plus près du corps 2 que du corps 1. Par conséquent, u est plus petit que 44, et, par 


conséquent, ( 1 — =) est positif. Le signe de la dérivée seconde correspond donc à celui de B, qui lui-même est 
Ua ui 
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déterminé par la valeur de i» : si cette quantité est supérieure à 1, alors B est négatif, alors que, dans le cas contraire, B est 
positif, ce qui implique que le point est instable. Le point de Lagrange L est situé au-delà du corps 2. La force totale 
(graitationnelle plus centrifuge) s'exerçant dans cette région est d'abord tournée vers le corps 2 quand on est proche de celui- 
ci, puis s'annule en L et est après dirigée à l'opposé de Lp. Au point tel que u» vaut 1, la composante de cette force, selon 
l'axe reliant les deux corps, est donnée (à une constante multiplicative positive près) par 

1 1 
ðf x —(1-q) |5 -u| -qg |3- ul 
w Ug 

avec, ici, 

Ua = l. 

LESI = d, 
soit 

T 

af xX g4 En q) 
Cette quantité étant strictement positive, le point u> = 1 de l'axe est situé au-delà du point L2. En conséquence, au point L>, U> 
est inférieur à 1, donc B est positif, donc le point est bien un point selle, ce qui assure de son instabilité. Une démonstration 
strictement analogue peut être faite pour le point L3, ce qui achève la démonstration de leur instabilité du fait de leur caractère 
de point selle. 
Temps caractéristiques en L et L, pour les systèmes à grande hétérogénéité de masse 

[modifier le code] 

Une des applications les plus importantes de l'instabilité des points de Lagrange, L4 et L>, réside dans le fait que des 
satellites artificiels peuvent être envoyés en ces points du système Terre-Soleil (voir ci-dessous). Pour de tels 
satellites, des corrections de trajectoires régulières doivent être appliquées afin de conserver le satellite au voisinage 
du point. Ce temps caractéristique peut être évalué dans le cas où le rapport de masse des deux corps du système 
est élevé. Dans ce cas, le temps caractéristique Va d'instabilité est donné par 

1 jy 

Y 2rÿ1+2V7 
où T est la période orbitale du système. Dans le cas du système Terre-Soleil, où T est légèrement supérieur à 365 
jours, le temps caractéristique d'instabilité est alors de 23 jours et 4 heures. 
Par ailleurs, la composante stable de la trajectoire se fait à la pulsation 

WTraj = wy2V7 — | 
soit, de façon équivalente, avec la période 

T 
Fis T 
V2V7—1 
ce qui, dans le même cas de figure que ci-dessus, donne une période de 176 jours. 
Démonstration 
L'équation donnant les \aleurs propres du système est toujours 
4 2112 2 — 

avec, pour les points L4 et Lp, 

Mszy — 0 

2 
() vx = —20 y ur Jw A 
qg l-q 

te= | = + — 1 | 
en se restreignant aux termes d'ordre le plus bas en g, u4 vaut 1, et up est déterminé par la relation donnée par le premier 
tableau de cette page. On a ainsi 

2 2 

Qy =W (3+1—-1) = 3" 
L'équation polynomiale devient alors 

À — Qu N — Qu = À 
dont les solutions sont 

X = (142V7) «À 
La solution positive à cette équation indique que les écarts au point d'équilibre croissent exponentiellement au cours du temps 
selon la relation 

5X x et 
avec 
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À =wy1+2V7=— TV +947 


Le temps caractéristique associé est donc 


~ 0,0635 T 


T 
N EEN 


soit, comme annoncé, un temps caractéristique de l'ordre de 23 jours pour les points de Lagrange de la Tere. 


Pr 
Î 


De la même façon, il existe des trajectoires périodiques dont la pulsation est donnée par les racines complexes de l'équation, 
soit 


WTraj = wy 2T — Í’ 


c'est-à-dire une période de 


{Traj m= D A = 0,4827 T 


2/7 — 1 


ce qui correspond à un temps de presque six mois pour les points de Lagrange de la Terre. 


La structure des orbites en présence d'instabilité [modifier le code] 


Une fois les valeurs propres d'un point instable connues, une trajectoire au voisinage d'un point de Lagrange va être 
une combinaison linéaire des vecteurs propres associés aux valeurs propres. En notant 4’ l'une de ces valeurs 
propres, le vecteur propre associé a pour composantes 


1 
V. = Hi 
Àifli 


m= _—_ (= SE à) 


5Y 
"ama 


où les quantités a; sont des nombres quelconques déterminés par la valeur des òX, ðY et de leur dérivée à un instant 
donné. Dans le cas des trois points de Lagrange instables, le déterminant de la matrice des dérivées secondes est 
négatif, ce qui implique que le discriminant de l'équation du second degré en A2 possède des racines réelles de 
signes opposés, et que, au terme, les valeurs propres recherchées sont deux nombres imaginaires purs opposés et 
deux nombres réels opposés. Une trajectoire générique comprend donc, dans le plan orbital, une composante 
périodique (liée aux racines imaginaires pures), une composante amortie (liée à la racine réelle positive), et une 
composante instable. Pour une position òX, ðY donnée, il est toujours possible de choisir une vitesse telle que les 
deux vecteurs propres aux racines réelles ne contribuent pas à la solution correspondante. La trajectoire obtenue est 
alors périodique, la période étant donnée par la racine complexe. Une telle solution n'est cependant pas stable. Un 
écart de trajectoire infime va en réalité rajouter à la trajectoire une composante instable, qui Va, peu à peu, éloigner 
la trajectoire de sa composante périodique. On dit que la trajectoire obtenue n'est pas dynamiquement stable. Ceci 
est une généralisation du fait qu'un objet situé exactement sur un point de Lagrange instable est dans une situation 
instable : un petit écart à cette position d'équilibre, inéluctablement généré par les perturbations causées par les 
autres corps du système, finira par éloigner l'objet de sa position initiale. La même chose se produit pour des 
trajectoires situées autour du point d'équilibre instable. 


Pertinence du concept [modifier le code] 


Le calcul ci-dessus fait référence à une configuration où les deux corps du système sont en orbite circulaire. 
Néanmoins, le concept de point de Lagrange prévaut pour tout type d'orbite, y compris elliptique. On peut donc 
définir ces points dans tout système à deux corps liés gravitationnellement. Par contre, les trajectoires, stables ou 
instables, autour des différents points de Lagrange dépendent explicitement de la circularité ou non de l'orbite des 
deux corps du système. 


Dans le système solaire 1mocifierie code] 


2% Article connexe : Liste d'objets du Système solaire situés à un point de Lagrange. 


Troyens [modifier le code] 


À Article détaillé : Troyen (astronomie). 


Les points L4 et L5 sont généralement stables, aussi on y trouve de nombreux corps naturels, dits troyens : 
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e dans le système Soleil-Jupiter, on recense (en 2011) environ 5000 astéroïdes aux points L; et Ls 

e dans le système Soleil-Neptune, huit 

e dans le système Soleil-Mars, quatre 

e dans le système Saturne-Téthys, les points L4 et L5 sont occupés par Télesto et Calypso, respectivement 
e dans le système Saturne-Dioné, Hélène et Pollux occupent ces points 


Curieusement, il semblerait que le système Soleil-Saturne ne soit pas en mesure d'accumuler des troyens du fait des 
perturbations joviennes. 


Dans le système Soleil-Terre, on connait (depuis peu) un troyen au point L,, l'astéroïde 2010 TK7, qui mesure 300 
mètres de diamètre. Certains astronomes soulignent que cet objet pourrait représenter un risque comparable aux 
géocroiseurs a Ces auteurs proposent également que l'impacteur supposément à l'origine de la Lune (Théia) aurait 
stationné un temps sur le point L4 ou Ls et accumulé de la masse avant d'en être éjecté sous l'action des autres 
planètes. 


Applications [modifier le code] 


Les points L4 et L sont des équilibres instables, ce qui les rend utilisables dans le cadre de missions spatiales : on 
n'y trouve pas de corps naturels, et un équilibre dynamique peut y être maintenu pour une consommation de 
carburant raisonnable (le champ gravitationnel étant faible dans leur voisinage). 


Les principaux avantages de ces positions, en comparaison des orbites terrestres, sont leur éloignement de la Terre 
et leur exposition au Soleil constante dans le temps. Le point L4 se prête particulièrement à l'observation du Soleil et 
du vent solaire. Ce point a été occupé pour la première fois en 1978 par le satellite ISEE-3, et est actuellement 
occupé par les satellites SoHO et Advanced Composition Explorer. Le point L est à l'inverse particulièrement 
intéressant pour les missions d'observation du cosmos, qui embarquent des instruments de grande sensibilité devant 
être détournés de la Terre et de la Lune, et fonctionnant à très basse température. Il est actuellement occupé par les 
satellites Herschel, Planck, WMAP et Gaia. 


Il a été un temps envisagé de placer un télescope spatial au point L4 ou Ls du système Terre-Lune, mais cette option 
a été abandonnée après que des nuages de poussière y ont été observés. 


En science-fiction [modifier le code] 


En science fiction, de par leur stabilité, les points L4 et L5 du système Terre-Lune abritent souvent de gigantesques 
colonies spatiales. Les auteurs de science-fiction et de bande dessinée aiment placer une Anti-Terre au point L3. 
Cette idée date d'avant la physique newtonienne, qui démontre qu'elle est tout à fait irréaliste. Le point de Lagrange 
n'a d'intérêt que pour un objet de masse négligeable par rapport aux deux éléments du système, ce qui n'est pas le 
cas d'une planète jumelle. Sans compter l'instabilité du point en question : une Anti-Terre placée en L3 tiendrait sur 
son orbite un peu moins d'un mois dans les meilleurs cas l'éf. souhaitée]. 


Notes et références jmodifierle code] 


1. 1 [Essai sur le Problème des Trois Corps htip:/www.Itas-vis.ulg.ac.be/cmsms/uploads/ÆFile/Lagrange essai _3corps.pdf ] 
2. T Géométrie de Roche, Jean-Marie Hameury, Observatoire de Strasbourg|1] 
3. 1 Bernard Bonnard, Ludovic Faubourg et Emmanuel Trélat, Mécanique céleste et contrôle des véhicules spatiaux, Berlin, 
Springer, coll. « Mathématiques et applications », 2005, XIV-276 p. (ISBN 978-3-540-28373-7, notice BnF n° FRBNF40153166), 
p. 73 (lire en ligne ) sur Google Livres (consulté le 25 juillet 2014) 
4. T Si l'on définit g comme le rapport de la plus petite masse à la totale, seules les valeurs de q inférieures à 0,5 ont un sens, 
puisque les valeurs plus grandes correspondent au rapport de la plus grande masse à la masse totale. 
5. 1 Do gravity holes harbour planetary assassins? 
e Grégory Archambeau, Étude de la dynamique autour des points de Lagrange, Orsay, Université Paris XI, 2008, X- 
114 p., lire en ligne [EDF] sur [2] 
e (en) Richard Fitzpatrick, Analytical Classical Dynamics, cours donné à l'Université du Texas à Austin Voir en ligne 
(versions pdf et html) 


Voir AUSSI Įmodifierle code] 


Articles connexes [modifier le code] 
| A | Sur les autres projets Wikimedia : 
e Liste d'objets du Système solaire situés à un point de Lagrange 


9 Point de Lagrange, sur Wikimedia 
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e Notices d'autorité : Bibliothèque nationale de France Librairie du zii point de Lagrange, sur le Wiktionnaire 
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v-:d'm Orbites 

Box (en) - Capture - Circulaire - Elliptique : Elliptique élevée : 

Évasion (en) - Orbite de rebut - Trajectoire hyperbolique (en) : Inclinée : 
Non inclinée (en) - Osculatrice : Trajectoire parabolique - Orbite 
d'attente : Synchrone (Semi-synchrone (en) - Sous-synchrone (en) ) 


Général 


Types d'orbite Géosynchrone : Géostationnaire + Héliosynchrone : Terrestre basse : 
Géocentrique Terrestre moyenne : Terrestre haute - de Molniya - Quasi équatoriale : 
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Orbite lunaire : Polaire + Toundra : Deux lignes NORAD 


Non Aréosynchrone (en) - Aréostationnaire (en) © Halo (en) - Lissajous - 
géocentrique Lunaire - Héliocentrique + Héliosynchrone 


3 Inclinaison - () Longitude du nœud ascendant : e Excentricité : ;,; 


Classiques Argument du périastre : ;4 Grand axe : M, Anomalie moyenne à 
n l'époque 
Paramètres , | | OTN | 
y Anomalie vraie : } Petit axe : g Excentricité -+ z? Anomalie 
Autres excentrique - L Longitude moyenne (en) ] Longitude vraie (en) : P 
Période orbitale 


Transfert bi-elliptique (en) : Bilan des modifications du vecteur vitesse (en) - Orbite de transfert 
géostationnaire - Assistance gravtationnelle + Trajectoire à incidence nulle (en) : Transfert de 
Manœuvres Hohmann : Transfert à faible énergie (en) - Effet d'Oberth (en) : Correction d'inclinaison (en) : 
Correction du phasage (en) - Rendez-vous - Transposition, docking, and extraction (en) : 
Évitement de collision (vaisseau spatial) (en) 
Apsides - Système de coordonnées célestes - Rétrograde - Époque - Éphéméride + Système 
de coordonnées équatoriales - Trace au sol : Interplanetary Transport Network : Lois de Kepler : 
Mécanique orbitale Point de Lagrange : Problème à N corps : Orbit equation (en) : Vitesse orbitale : Vecteurs 
d'état d'orbite (en) : Perturbation - Specific orbital energy (en) - Specific relative angular 
momentum (en) 


Liste d'orbites 


€ Portail de l'astronomie EN Portail de l’astronautique 


Catégories : Mécanique céleste | Joseph-Louis Lagrange 
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